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Abstract 


Carbon-dioxide (CO2) is the main contributor to anthropogenic global warming, and the timing of 
its peak concentration in the atmosphere is likely to govern the timing of maximum radiative forcing. 
It is well-known that dynamics of atmospheric CO2 is governed by multiple time-constants, and here 
we approximate the solutions to a linear model of atmospheric CO2 dynamics with four time-constants 
to identify factors governing the time-delay between peaks in CO2 emissions and concentrations, and 
therefore the timing of the concentration peak. The main factor affecting this time-delay is the ratio 
of the rate of change of emissions during its increasing and decreasing phases. If this ratio is large 
in magnitude then the time-delay between peak emissions and concentrations is large. Therefore it is 
important to limit the magnitude of this ratio through mitigation, in order to achieve an early peak in 
CO2 concentrations. This can be achieved with an early global emissions peak, combined with rapid 
decarbonization of economic activity, because the delay between peak emissions and concentrations is 
affected by the time-scale with which decarbonization occurs. Of course, for limiting the magnitude of 
peak concentrations it is also important to limit the magnitude of emissions throughout its trajectory, but 
that aspect has been studied elsewhere and is not examined here. The carbon cycle parameters affecting 
the timing of the concentration peak are primarily the long multi-century time-constant of atmospheric 
CO2, and the ratio of contributions to the impulse response function of atmospheric CO2 from the infinite 
time-constant and the long time-constant respectively. Reducing uncertainties in these parameters can 
reduce uncertainty in forecasts of the radiative forcing peak. 
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1 Introduction 


As countries agree on commitments towards a new international climate treaty to be decided in 2015 


forcers 


(Forster et al 


iUNFCCC (2014a b )), the se will include mitigation of not only carbon-dioxide (CO 2 ) but also other climate 


(I2OO7II : 


UNFCCC ( 2014c)). CO 2 is , and is likely to remain, the largest contribution to radiative forcing 


Mvhre et al. 


(1201311 1. Limiting long-term warming requires limiting the growth in 


global CO 2 emissions, and eventually reducing these emissions. If the present increasing trend in global 
CO 2 emissions is eventually reversed so that an emissions pe ak occurs, th e corre s ponding peak in concentra - 


tion will be delayed because of its long atmospheric lifetime \Allen et al 


(I2OO9II : 


Meinshausen et al 


(I2OO9II : 
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Mianone et al 


(l2008f) l. A CO 2 concentration peak would be a significant event for global climate: it would 


govern the maximum contribution of CO 2 emissions to radiative forcing. Furthermore, assuming that CO 2 
continues to be the major contribution to radiative forcing, then its peak concentration will strongly influence 
the magnitude and timing of peak global warming. 

The Earth’s C O 2 cycle is complex, i nvolving multiple reservoirs that maintain exchanges occurring at very 


different rates ( Archer et al 


(Il997h : 


Cox et al. 


( 200 ( 111 : 


Falkowski et al. 


excess CO 2 is by the surface ocean and land biosphere \Pierrehumhern (|2014l ll. Progressively slower pro¬ 


( 200 c 11. The most rapid uptake of 


cesses involve mixing with the deep-ocean, reduction of ocean acidity due to dissoluti on of carbonates, an d 


uptake of excess atmosphe r ic CO 2 via reaction with CaCOs or silicate rocks on land \Archer et al. 


Archer and Brovkin (1200811 : 


Archer et al. 


(I1997II : 


( 200911 1. The last two processes require many tens of thousands of 


years so that, on the timescales of the next few centuries, their contributions can be effectively neglected. 
Equivalently their effects can be treated as occurring with infinite time-constant. Accurate characterization 
of the different processes involved, in order to describe the fate of excess atmospheric CO 2 , requires coupled 
climate-carbon-cycle or Earth-system models; such models ha v e bee n employed to describe effec ts of mit¬ 


igation scenarios on CO 2 in the atmosphere \Petoukhov et al. 


( 200511 : 


Friedlinastein et al. 


( 200611 1. As the 


mitigation of CO 2 emissions unfolds, these and similar models will play important roles in estimating the 
consequences for atmospheric CO 2 , including the timing and magnitude of its peak concentration. 


This paper solves a linear model of atmospheric CO 2 with four time-constants \Joos et al. 


( 201311 1 to under¬ 


stand the factors controlling the time-delay between peaks in emission and concentration and therefore the 
timing of the concentration peak. Previous studies have described th e relationship between m i tigation and 


warm i ng, and highlighted 


(120091); 


Allen and Stocker ( 2014 1; 


the importance of rapid mitigati on (for e.g. \Socolow and Larrn (120071) ; 


Allen et al 


Huntinaford et al. 


(I2OI2I II. Here we focus specifically on solving the model 


of atmospheric CO 2 analytically, to identify some of the important factors controlling the time to the con¬ 
centration peak of C 02 . 


2 Models of emissions and carbon cycle 


2.1 Carbon cycle model 


Joos et al. 


(I2OI3I ) estimated the impulse response of CO 2 in the Earth’s atmosphere, by averaging a group 
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of Earth System Models. They estimated a mean response with four time-constants 


h (t) = 0.276e"*/'‘-3° 0.282e“‘/^®-^ -f 0.224e"*/3®'‘ -f 0.217 


( 1 ) 


which we apply, and compute atmospheric concentration using 


{t) = / m{z)h{t — z) dz + Up 

Jq 


( 2 ) 


where m (z) is anthropogenic emissions in concentration units, i.e. the rate of increase in concentration in 
the hypothetical case of infinite atmospheric lifetime. The constant term up is preindustrial concentration. 
Concentration is noted in parts per million (ppm). Atmospheric emission of CO 2 in the year 2013 was 
36 X 10^^ kg, equivalent to 4.5 ppm^ We furthermore describe the impulse response function symbolically 
as 


K 


h{t) = Y^ 


IjLie 


-t/T 


( 3 ) 


with = 1 and {ti,T 2 ,T 3 ,T 4 } = {4.4, 36.5, 394, 00 } following the results of 


Joosetal 


(l2013h . The 


infinite time-constant approximates the effects of very slow processes involving buffering of ocean acidity by 
dissolution of carbonates and the uptake of CO 2 in the weathering of rocks. 

There is un c ertain ty in the function h (t ), with different earth system models likely to yield different results 


{Joosetal 


(20131)). While the present paper does not characterize this uncertainty. Section 3.3 considers 


the influence of changing those parameters in this 4-time-constant model that are shown to affect the timing 
of the concentration peak. 

2.2 Emissions model 


The model of emissions m (t) is very simple, and chosen so as to describe emissions using a few different pa- 
rameters that can be readily interpreted. The emissions model is m (t) = mo (1 -I- g-{t-to)/T^ 


^Seshadrn (l2015al la)). with mo being present emissions, r the growth rate of gross global product (GGP), and 


to denoting the present time. It is assumed, for the purposes of the emissions model used here, that GGP 
increases for tg years from the present at constant rate r, after which it remains fixed. The term e-d-to)lrm 
describes the effect of decrease in emissions intensity of economic output, with ^ 00 corresponding to 
the absence of any mitigation, and smaller values of Tm. co r respon ding to rapid mitigation. This model of 


emissions has been borrowed from previous work ( Seshadri ( 2015 JLd)- 


CO 2 had infinite lifetime, emissions in the year 2013 would have increased atmospheric concentration by 4.5 ppm. 
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3 Results 


3.1 Peak atmospheric concentration of C02 


With emissions m (t) the atmospheric concentration is u{t) = up + m (z) h(t — z) dz, and differentiating 
this equation we obtain for the rate of change of concentration that u' {t) = m{t) + f* m {z) dz, 

where we have used the fact that h (0) = 1 from equation ([T]). Using equation ([2]) we obtain that u' (t) = 
m (t) — /p (z) dz. The last integral can be evaluated by parts to give (z) dz = 

Tim {t) — Ti f* (z) dz, and substituting this yields 


K ft 

' (t) = ^ (z) dz 

Jo 


(4) 


for the rate of change of concentration. This result has been derived previously in 


Seshadri (2015c). 


For short time-constants ti = 4.4 years and T 2 = 36.5 years for which t ^ ti, i = 1,2, we can approximate 
the integral (z) dz by {t) . For the very long time-constant that is represented by r 4 = oo 

the integral becomes /p e^/'^'m' (z) dz = /q m' (z) dz = m {t) . Hence the rate of change of concentration 
becomes approximately 


u'(t) = {fiiTi + ^j. 2 T 2 ) m'(t) + Use {z) dz + 

Jo 


(5) 


We denote the time to the emissions peak as ti, the time to the concentration peak as t 2 , and the time-delay 
as St = t 2 — ti. Then the integral in equation ([S]) above can be described as the sum of integrals from 0 to 
ti, and from ti to 12- Furthermore, defining weighted rates of change of emissions 


and 


(z) dz 


Iti 


we obtain the formula for the rate of change of concentration at time t = t 2 


( 6 ) 

(7) 


u' {t 2 ) = (miTi -I- ^. 2 T 2 ) m' (t 2 ) + M3T-3m^„ye - 1^ +/i3'r3TOjj„_^e -|-^4TO (^ 2 ) 

( 8 ) 


Before proceeding we identify the factors influencing whether the concentration will reach a peak value and 
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eventually decline, as opposed to continuing to increase to an asymptotic maximum. Consider emissions 


scenarios where the emissions peaks and then declines to zero so that, eventually. 
In that case only the large but finite time-constant T 3 = 394 years in the model of 


m' (t) = 0 a.nd m (t) = 0 . 


JoosetaL 


( 2013 1 plays a 


role in the sign of u' (t). The value of this quantity then becomes, for t > ti, approximately 




(9) 


which is negative if 


, , - 1 
“ '^av,d > ''^av,i^t/T3 _ gti/xa 


( 10 ) 


requiring the rate of decrease of emissions to be sufficiently large in magnitude. The larger the average rate 
of increase in emissions, and the longer the increase persists, the more stringent is the condition on the 
subsequent decrease in order for concentrations to eventually decrease. Therefore an early peak in global 
emissions can help stabilize atmospheric CO 2 concentrations, as would be expected. 


We now seek an expression for the time t 2 to the CO 2 concentration peak, assuming that mitigation occurs 
rapidly enough for one to occur. The condition for such a peak is that the rate of change of concentration 
vanishes. Hence, from equation ([5]) above 


u' (t 2 ) = (Min + M 2 T 2 ) m' {t2)+H3T3m'^v,i(^ - l^+fi3T3m'^^ j^e (^ 2 ) = 0 

( 11 ) 

and, writing m (12) = rn'^ti + (^2 — ^i); where m' and are average rates of change during emission’s 

increasing and decreasing phases respectively, and collecting terms involving t 2 yields 


Mane - m'^v,d) +^J-Am'^t 2 = -M 3 T 3 TOL.d-M 4 (m' - mj^) ti-(Min -k M 2 T 2 ) m' (^ 2 ) 

( 12 ) 

We can neglect the last term on the right side of equation since m' { 12 ) is comparable in magnitude to 
d; Eis will be shown, but Min + M 2 n Man- Hence the time ^2 to the concentration peak in the model 
is governed by approximate equality 

Mane"*"/""" - m'^v,d) - M4TOrft2 = Manw).^^^ -k M4 (m' - to^) ti (13) 

The approximate equality in equation (11311 has been verified in Figure Ic, for emissions scenarios shown 
in Figure la and corresponding concentration graphs plotted in Figure lb. It is therefore confirmed that 
the short time-constants can be approximately neglected while studying the concentration peak. The above 
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equation has to be solved numerically because neither the exponential nor the linear term in t 2 on the left 
side of the expression can be neglected. That this is the case is shown in Figure Id, which plots these two 
terms. Neither term is dominant. 


However in order to understand qualitatively the factors influencing the time t 2 , let us imagine that t 2 is so 
large that the second term in the left side of equation m, —H 4 m'^t 2 , were dominant. Then the solution 
would be given by 


t 2 = T3 


-m' 


+ 1 


^ _l_ M3 
T3 fj,4 -m'^ 


(14) 


The time to the concentration peak increases with the rate of increase of emissions during their growing phase. 
It decreases with the rate of decrease of emissions during their declining phase. The time to the concentration 
peak increases with the ratio ^ 44 /^ 3 , describing the ratio of the impulse response of atmospheric CO 2 coming 
from the infinite time-constant T 4 and finite but long time-constant T 3 respectively. It increases with the 
time ti to the emissions peak, and in fact the influence of ti can be described in terms of dimensionless 
ratio ti/ra. However the left side of equation (fOll increases with H, because of the exponential term, so the 
influence of H is not as strong as equation m, which neglects the influence of this term, would suggest. 
Similarly the influence of T 3 is not as strong as equation (fHll would suggest because the exponential term in 
equation ca also increases with this quantity. 


We can write equation m in terms of dimensionless variables il'^'av d) parameter ^ 4 /^, 3 , by 

writing it in equivalent form 


Tae 


-tllT-i 


\m'av,d 



M4 TOd ^ 

- } -C 2 

M3 


= T3 


M 3 


TO,- 


^av,d 




(15) 


The above model depends also on variables to'/to(j„^ and Recall that these are ratios of un¬ 
weighted rates of change and weighted rates of change, i.e. weighted by It is shown later that ^ 

is approximately constant across scenarios (also see Figure 3). 

Figure 2 plots the time-delay between peak emissions and concentrations versus the dimensionless variable 
'^av il'^av d^ the emissions scenarios plotted in Figure 1. The time-delay increases with the absolute 
value of this ratio. Furthermore, for scenarios with relatively short e-folding mitigation timescale, less than 
about 40 years, the relationship is approximately linear for the family of scenarios considered here. 
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Figure 1: Verification of approximation in equation (1131) and demonstration that neither term on the left 
side of this equation is dominant for a wide range of CO 2 emissions scenarios: (a) emissions pathways; (b) 
corresponding graphs of CO 2 concentration; (c) each side of equation (13), with left side in abscissa and right 
side in ordinate. The dashed line shows where abscissa and ordinate are equal; (d) first and second terms on 
the left side of equation dni), /i = and /2 = , showing 

that neither term is dominant. Dashed line shows were abscissa and ordinate are of equal magnitude. Colors 
in panels correspond to different e-folding mitigation timescales in the respective emissions scenarios (see 
legend in Figure 2). 



































Figure 2: Time-delay 5t between peak emissions and concentrations versus the ratio of the weighted-averaged 
rate of increase and decrease of emissions, defined in equations ([5]) and © respectively. The time delay 
increases with the absolute value of this ratio. Note the reversed axis in the abscissa. For short e-folding 
mitigation timescales, the relationship is approximately linear. Results are shown for the scenarios in Figure 
1, and colors in the plot correspond to the colors in Figures la-d. 

3.2 Influence of mitigation parameters, and importance of e-folding mitigation 
timescale 


Here we consider the effects of parameters of our specific emissions model on the rates j, ^ , m' and 
, and thereby on the solution to equation (uni). This discussion is of more general interest beyond this 
particular model because the parameters - the GGP growth rate, the e-folding mitigation timescale at which 
emissions intensity decreases, and the time to stabilization of GGP - can be easily interpreted. Gonsidering 
first the rate of change of emissions during its increasing phase 

, = ’HM. (16) 

H'dz f. 

Approximating (1 -|- r)* = e’'* the formula for emissions becomes m (G) = If r > I/t™ 

then ti = tg, i.e. emissions peaks when GGP is maximum. Otherwise ti = to, the present. In the first case 
m{ti) = and m' = whereas in the second case m(ti) = toq and 

TO- = mo/to- 

In the following we consider only the case where r > I/t^ because global emissions of CO 2 are expected to 
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continue to increase for a while. For its decreasing phase 


, m{t 2 )-m{ti) 

ruj = — 


i::dz 


t 2 — ti 


( 17 ) 


Then m (^ 2 ) = m (<i) e~ so that m'^ = —m (<i) ( 1 — e~ ) / (^2 — ^i)) which becomes — 


For the weighted rate of increase of emissions between 0 and <1 , we decompose the integral in the numerator 
of equation dS]) into that between 0 and to and between to and ti. Using m' (t) = 
between to and ti we obtain 


/ ^-^ / {r ^ 1 f 

Tm T 3 ^ 


9 .^3 _ g t3 


) 


T3 (e*9/-3 - 1) 

where (z) dz. Likewise using m! (t) = —between ti and t 2 


moe 






(t3 - Tm) (e*2/'^3 - 


(18) 


(19) 


Let us compare the rates mj; and ^ by considering ratio 


e^3 


^av,d 


_ *2 ~*g I *2 
1 — e Tm ^ T 3 


t 2 - ta 


1 — e 


^ t 2 —*g 

(t3 - Tm) 6^3 ( e -3 


- 1 


( 20 ) 


which, assuming Ijrm^ I/T 3 because the mitigation timescale is generally much shorter than this time- 
constant of CO 2 , simplifies to 

( 21 ) 


* 2 -tg 

(t 3 - Tm) e -3 - 1 


In case ° <C 1 the exponential term in equation m simplifies to 1 -|- ’’ , in which case the ratio 


^2 —t 
T3 

approximates to 




;,d ^ "^3 


m' ra - 


( 22 ) 


which, for cases where T 3 ^ Tm, approximates to 1. Hence the value of weighted average while not 

exactly equal that of the average rate of change of emissions in its decreasing phase, closely approximates 
the latter, especially in cases of rapid mitigation when, in addition, the time-delay between peak emissions 


to) 
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Figure 3: Plots of weighted (by versus unweighted rate of change of emissions: (a) increasing phase; 

(b) decreasing phase. Dashed lines show where abscissa and ordinate are equal. An approximate linear 
relationship holds between the two sets of variables. 


and concentrations is short compared to time-constant A near-equality between these two variables, 
weighted and unweighted, is seen in Figure 3b. Similarly Figure 3a shows the relationship between weighted 
and unweighted rates of change for the phase when emissions are increasing. While departure from equality 
is larger, the relationship is still approximately linear. 

Although these weighted and unweighted rates are not the same, because of their similarities we can examine 
ratio , which is more tractable, to understand qualitatively what controls the behavior of ratio 

The former ratio can be written as 


m. 


t2 - tn 


1 — e 


(23) 


A longer time t 2 to the concentration peak by itself can increase the above ratio, because the graph of 
emissions is convex in its decreasing phase, so that its slope is decreasing. The influence of time tg to 
peak emissions is weak because of two countervailing influences: shorter tg increases the magnitude of the 
numerator as well as that of the denominator above. The main influence on this ratio is that of mitigation 
timescale . Short decreases the magnitude of this ratio. 

The strong influence of the mitigation timescale on the ratio ^ is seen in Figure 4. Short mitigation 

timescale, corresponding to rapid mitigation of emissions intensity, is therefore essential to limit the absolute 
magnitude of this ratio and therefore assure a short time-delay. 

Note that the ratio in equation (I23|l does not depend on the GGP’s growth rate. While this rate influences 
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Figure 4: Ratio versus e-folding mitigation timescale. Note the reversed axis on the ordinate. 

Short mitigation timescale decreases the absolute value of this ratio. 

peak emissions of CO 2 , it affects the average growth rate and decrease of emissions in the same manner 
and hence is not a factor in this ratio and consequently in the time-delay between peak emissions and 
concentrations. 


3.3 Carbon cycle uncertainties 


As indicated earlier the carbon-cycle parameters affecting the time to the concentration peak are the multi¬ 
century time-constant T 3 and the ratio /r 4 //r 3 , describing the ratio of the impulse response of atmospheric CO 2 
from the infinite time-constant T 4 and the long time-constant r 3 respectively. Figure 5 plots the influence of 
these parameters on the time-delay. These are the uncertainties in the carbon cycle that must be constrained 
in order to constrain forecasts of the timing of peak concentrations of CO 2 . 



4 Conclusions and Discussion 


T he results presente d here are based on approximating the linear carbon cycle model with four time-constants 


of 


Joos et al. 


(12013), wi t h one t ime-constant being infinite because a fraction of emitted CO 2 persists for a 


very long time iArchen (|2005l l: \Archer and Brovkin (1200^). We identified the main factors governing the 
time-delay between peak CO 2 emissions and concentrations. 
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Figure 5: The main carbon cycle parameters affecting the time-delay between peak CO 2 emissions and 
concentrations: (a) relation between time-delay and ^ 4/^3 describing the ratio of the impulse response of 
atmospheric CO 2 from the infinite time-constant T 4 and the long time-constant T 3 respectively; (b) relation 
between time-delay and long time-constant T 3 . In each panel, different curves correspond to different fixed 
values of the other parameter. The emissions scenario used has tg — to = 60 years, r = 0.015 %/year, and 
Tm = 50 years. 


On the emissions front, the main factor is the e-folding timescale with which the mitigation of emissions 
intensity of GGP ("decarbo nization") occurs. This can be viewed as the inverse of the corresponding mitiga¬ 


tion rate iSeshadri 


(l2015al lb Short decarbonization timescale leads to short time-delay between emissions 


and concentration peaks. Therefore achieving decarbonization rapidly is important to achieving an early 
peak in CO 2 concentrations. 


The time-delay between peak emissions and concentrations is not sensitive to the time to peak emissions. 
However an early emissions peak will facilitate an early concentration peak. 

The growth rate of economic output is an important factor in peak emissions. However as discussed here 
it does not affect the time delay between peak emissions and concentrations, because it has the same effect 
on the rate of increase of emissions and the rate of decrease of emissions, whose ratio governs the time- 
delay. Therefore in this model where peak emissions corresponds to where GGP stabilizes, the timing of 
peak concentrations is not affected by this growth rate. Nevertheless it is obviously important in consid¬ 
ering the magnitude of peak emissions, and faster economic growth has to be accompanied by more rapid 
decarbonization. 


The important non-dimensional parameter in our discussion has been the ratio of the rate of increase of 
emissions and the rate of decrease of emissions. Limiting the magnitude of this ratio will help achieve an 
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early peak, and this can be accomplished by keeping the mitigation timescale short. 

With respect to the atmospheric cycle of CO 2 , the influential parameters are the long but finite time-constant 
T 3 that occurs on century-scales, and the factor ha /^3 describing the ratio of the impulse response function 
of atmospheric CO 2 from the infinite time-constant T 4 and the long time-constant T 3 respectively. The time- 
delay increases with these parameters. The short time constants ri and T 2 occurring on decadal scales or 
less play a small role in the long-term dynamics of atmospheric CO 2 , and uncertainties in their values are 
correspondingly less important for forecasting the concentration peak. 

In summary it is important to constrain these carbon cycle parameters, in addition to achieving an early 
mitigation peak, as well as implementing decarbonization on short timescales. 
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